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The effect of mutual interaction between dark solitons and dispersive waves is investigated nu¬ 
merically and analytically. The condition of the resonant scattering of dispersive waves on dark 
solitons is derived and compared against the results of numerical simulations. It is shown that 
the interaction with intense dispersive waves affects the dynamics of the soltons strongly changing 
their frequencies and accelerating or decelerating the solitons. It is also demonstrated that two 
dark solitons can form a cavity for dispersive weaves bouncing between the two dark solitons. The 
differences of the resonant scattering of the dispersive waves on the dark and bright solitons are 
discussed. In particular we demonstrate that two dark solitons and dispersive wave bouncing in 
between them create solitonic cavity with convex “mirrors” unlike the concave “mirror” in case of 
the bright solitons. 


INTRODUCTION 

Optical dark solitons were predicted [T] and experi¬ 
mentally realized [2] several decades ago and since then 
attracted considerable attention |3H9] in optical commu¬ 
nity. The interest to the dark solitons is motivated by 
their fascinating fundamental properties and by possible 
applications for telecommunications 011]. Dark solitons 
exist in the normal dispersion regime and they are char¬ 
acterized by the dip on top of the flat intensity back¬ 
ground. Vectorial dark solitons uni as well as coupled 
bright and dark solitons were considered iniiE]. Also 
the dynamics of dark solitons in presence of higher order 
dispersive effects was well covered 0I1II1HTSI. Several 
studies have also addressed emission of dispersive waves 
(DWs) by dark solitons [13] [15] as well as emission of 
radiation by shock waves [16] propagating in normal dis¬ 
persion region. 

Discussing the interaction of dispersive waves and 
bright solitons it is relevant to remind that extensive 
studies were performed on scattering of DWs on bright 
solitons, with solitons properties remaining invariant [IS- 
[22]. In particular it was also demonstrated that strong 
resonant scattering of DWs on solitons can accelerate or 
decelerate the solitons with consequent frequency up- or 
downshifts I23l[2l|. This mechanism can be viewed as 
a method of all-optical switching [24ll26] or alternative 
technique for generation of broad and coherent super¬ 
continuum [25]. A very interesting scenario occurs when 
DWs emitted by one soliton interacts with a neighboring 
solitons through the dispersive waves trapped between 
the solitons. The multiple scattering of DW on solitons 
results in the mutual attraction of the solitons and even¬ 
tually to their collision m- Furthermore, these collisions 
can cause fusion of two solitons [28]. This phenomenon 
was proven to be responsible for the appearance of multi¬ 


ple soliton knot patterns [23l [29j |30] during the complex 
super continuum generation process. Very recently the ef¬ 
fect of strong variations of the solitons trajectories under 
the action of the DW was experimentally demonstrated 
[26] . An experimental observation of the solitonic cavities 
confining DW was very recently reported in m- 

This work aims to address the interaction of dark soli¬ 
tons with external DWs and to explore the possibility 
of manipulation of dark solitons by DWs. We will also 
demonstrate a possibility to create solitonic cavities made 
from dark solitons and DWs oscillating in between them. 

INTERACTION OF WEAK DISPERSIVE WAVES 
WITH DARK SOLITONS 

For a weak DW interacting with the bright soliton it 
is possible to construct an analytic theory predicting the 
location of the resonance frequencies [TH] [19] . Here we 
develop such a theoretical approach based on the per¬ 
turbation theory for the case of dark solitons interacting 
with weak DWs. 

Dynamics of the dimensionless amplitude u{z,t) in 
normal dispersion regime with significant third-order 
dispersion can be described by a generalized nonlinear 
Schrodinger equation 

idzu - f + \u\'^u = 0 (1) 

We seek for a solution in the form of a dark soliton per¬ 
turbed by a sum of incident DW and scattered radiation, 
generated by interaction between DW and the soliton 

u{z,t) = + 'tpinc{z,t) + i;sc{z,t)] 6*"°^ 

where 2psoi = tanh(i4ot). Assuming that both DW 
and scattered radiation are small compared to the soliton 
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background Uq and retaining only the linear terms we can 
write the following equations for and 

1 % 

idzlpinc - V’inc - Q/Sadt-tpinc + ul {-Ipinc + V’t*„c) = 0 

( 2 ) 

1 i 

idz-ipsc - - -^l^adtP’sc + (2^/'L “ ul)^psc + V’LV’sc 

= - ‘^i'lps - «o)V’mc - (V’s - ulMnc 

(3) 

Right hand’s side of equation includes three driving 
terms: first one generates Cherenkov radiation of the soli- 
ton (which has been already analyzed in [l3]), while the 
other two correspond to FWM process between DW and 
the soliton. Seeking for a plane wave solution for (U) by 
setting 

-Ipinc = 


we arrive at the dispersion relation 

HcOinc) = ^Pzl^inc ± ‘^inc]J ( 4 ) 

The condition of the resonance excitation of the delocal¬ 
ized eigenmodes of the medium can be obtained by the 
analysis of the asymptotic of the eigenfunctions when 
where approaches Uq. The equation then turns into 
a version of (© with additional forcing terms and for the 
FWM resonance frequencies ujsc we can write 


^ ^sc 






( 5 ) 

To confirm the analysis above we have conducted a se¬ 
ries of numerical simulations. § was solved using split- 
step Fourier method with the initial condition 


u{0^t) = uot8inh.{uot)-\-Aex.p{—{t—ti)^/r^)ex.p{—iujinct) 


where A is the amplitude, ti is the initial position of the 
center and r is the width of DW. We started with a rel¬ 
atively weak dispersive wave with amplitude of A = 0.1 
reflecting off the soliton with = 1. A typical result of 
simulation is presented on Fig. As it is shown on Fig. 
[^b) around the collision point {z ~ 35) the frequency 
of the dispersive radiation changes. The new frequency 
of the DW is very well predicted by the resonance con¬ 
dition see Fig. [^c) showing the graphic solution of 
the resonance condition. To further test we have var¬ 
ied the frequency ujinc of the incident pulse of DW and 
measuring the scattering frequency uJsc at the end of sim¬ 
ulation distance. The numerically measured frequencies 
are shown in Fig. by red squares alongside with the 
prediction of the resonance condition shown by the black 
curve. As one can see, the resonance condition predicts 
the frequency of the scattered DW oogc very well. 


It should be noted here that the resonance condition 
predicts two resonant frequencies even if the dispersion 
charactersitics of the DW has only one branch. In our 
case the dispersion characteristics of the small excitation 
on the background of the dark soliton has two branches 
and the resonant scattering between the DW belonging to 
different branches of the the dispersion characteristics is 
also possible. However only one of the resonant scattering 
channels is efficient and in this paper we discuss only this 
kind of scattering. 



FIG. I: A relatively weak DW {A = O.IO) scattering on a 
dark soliton. (a) is the field intensity minus the background 
||r^( 2 :,t)|^ ~ (b) is the spectral density \u{z,(ju) \ in the rel¬ 

evant region of frequency domain, (c) displays a graphical 
solution to 1^, where the blue curve and the red line are, re¬ 
spectively, the left and the right hand’s sides of the equation. 



FIG. 2: Scattered frequency Usc as a function of incident fre¬ 
quency LUinc- Solid line is the theoretical prediction according 
to markers are found by numerical modeling. 
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FIG. 3: An intensive DW (A = 0.20) scattering on a dark 
soliton. (a) is the field intensity minus the background 
~ '^ol^ (b) is the spectral density \u{z,(ju)\ in the rel¬ 
evant region of frequency domain. 


FIG. 4: A very intensive DW (A = 0.30) colliding with a dark 
soliton triggering its collapse, (a) is the field intensity minus 
the background | t)|^ ~'^ol^ (b) is the spectral density 

\u{z,uj)\ in the relevant region of frequency domain. 


INTERACTION OF STRONG DISPERSIVE 

WAVES WITH DARK SOLITONS AND DARK 
SOLITONIC CAVITIES 

The interaction of an intensive DW with a soliton is a 
more complicated case with richer dynamics, which, un¬ 
fortunately, cannot be adequately described by the ana¬ 
lytics developed in the previous section. The main draw¬ 
back of the simple linearized theory is that we have as¬ 
sumed that the soliton does not change while interacting 
with DW, which is, strictly speaking, not entirely true 
even for a weak DW. Collision with more intensive DWs 
seems to have a substantial effect on the solitons and, 
as a consequence, on the spectral structure of the DW 
itself. For example, twofold increase in the amplitude of 
DW {A = 0.20) leads to DW pushing the soliton from 
its original trajectory as it is demonstrated on Fig. 
Importantly, the steering of the trajectory in case of the 
dark soliton takes place in the direction opposite to that 
of the incident DW. This is completely different from the 
case of the bright solitons that is steered toward the inci¬ 
dent DW. Looking at the process in the frequency domain 
(see Fig. |^b)) one can notice that scattered wave is con¬ 
siderably broadened spectrally. This can be explain by 
the fact that the soliton frequency is changing under the 
action of the incident DW and therefore different parts 
of the incident pulse scatter into different frequencies en¬ 
riching the spectrum of the DW. 

Further increase of the intensity of DW leads to a de¬ 
composition of the black soliton that manifest itself by 
a frequency continium in the spectral domain. Fig. 
demonstrates an outcome of the simulation in which an 
intensive DW with the amplitude A = 0.30 collides with 
the soliton causing disintegration of the soliton. 

A very interesting configuration arises when we launch 
a relatively weak DW in a free space bound by two dark 


solitons. Initial condition in this case is changed to 

u{0,t) = Ucvt{t) + Aex.p{-t^/T‘^)exp(-iLOinct) 

—tanh(i4o(t + ti)), t<0 

+1^0 tanh(i4o(t — ti)), t>0 

where ti is the soliton launch position. In this case two 
solitons act as boundaries trapping DW inside a soliton 
cavity. DW periodically bounces off the walls of the cav¬ 
ity while changing its carrier frequency in a periodical 
fashion, as it is illustrated in Figs. Since we are 

dealing with multiple scatterings, soliton trajectory can 
be significantly perturbed even by a relatively weak DW 
with amplitude A = 0.1 (see Fig. [^. There is a major 
difference from the case of two bright solitons creating 
a solitonic cavity nil ns EB. The two dark solitons 
create a “convex mirror” like cavity when reflecting the 
DW and not “concave mirror” cavity as produced by the 
bright solitons. 

CONCLUSIONS 

The resonant scattering of weak and strong dispersive 
waves on dark solitons is studied analytically and numer¬ 
ically. An anlytical formula predicting the frequency of 
the scattered wave was derived for the case of weak dis¬ 
persive waves. Numerical simulations show that the res¬ 
onant condition predicts the frequency of the scattered 
waves with very high accuracy. The interaction between 
the solitons and relatively intense DW was studied by 
means of pure numerical simulations. It is demonstrated 
that strong dispersive waves can affect dark soliton’s tra¬ 
jectory strongly which opens a way to control dark soli¬ 
tons by dispersive waves. An important observation is 
that the bright soliton trajectories bends against the di- 
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FIG. 5: A very weak DW {A = 0.05) trapped in a two soli- 
ton cavity, (a) is the field intensity minus the background 
~ '^o|? (b) is the spectral density \u{z,(ju)\ in the rel¬ 
evant region of frequency domain. 
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FIG. 6: A weak DW (A = 0.1) trapped in a two soli- 
ton cavity, (a) is the held intensity minus the background 
||i^(2:,t)|^ ~ '^o|? (b) is the spectral density \u(z,uj)\ in the rel¬ 
evant region of frequency domain. 


rection of the incident DW and the dark solitons trajecto¬ 
ries bend in the direction of the incident DW. It was also 
demonstrated that two dark solitons allow creation of a 
solitonic cavity being able to trap dispersive waves mak¬ 
ing them bouncing between the solitons behaving like a 
couple of convex mirrors. The radiation trapped between 
the solions results in the dispersive-waves mediated repul¬ 
sion of the solitons at the distances much larger than the 
characteristic size of the solitons. 
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